Method for Transmission Line Analysis

ABSTRACT

The present invention concerns methods for estimating one or more transmission properties of a telecommunications transmission line. Estimates of the line input impedance, Z in  at a frequency f, and the line capacitance Ĉ are made. An estimate is then made in dependence of Ĉ, the frequency f and a value Γ, where Γ satisfies the relationship Γ coth (Γ)=Z in ·jωĈ, and ω=2πf. Accuracy may be improved by least squares fitting a curve to a set of intermediate values and then generating an estimate from the coefficients of the curve.

TECHNICAL FIELD OF THE INVENTION

The present invention relates to the field of transmission line analysis.

BACKGROUND

It is of great interest for the operation of telecommunication networks to measure the properties of the transmission lines of the network. The results of such measurements can be used for e.g. detecting faults, locating faults, predicting certain faults which are about to occur, and estimating the suitability and capacity of a line for certain services, such as DSL.

A variety of methods and devices can be used to estimate line properties.

In patent application WO2005SE01619 (publication no WO2007050001) and the patent applications it refers to, a method is described wherein a signal containing certain frequencies may be sent through a line board, a resulting signal received and transmission properties estimated from the relationship between the signals.

The accuracy of estimates of line properties are often not sufficient for all purposes. Hence there is always a desire for more accurate, yet reasonably computationally efficient methods to make estimates of line properties.

SUMMARY

The present invention provides a new method to estimate one or more transmission properties of a telecommunications transmission line.

According to the method, estimates of the transmission line capacitance Ĉ, and the complex impedance Z_(in) of the transmission line at a frequency f are made. One or more transmission line properties are then estimated in dependence of Ĉ, f, and a value Γ that satisfies Γ coth(Γ)=Z_(in)·jωĈ where ω=2πf.

This has the advantage that better accuracy can be achieved than previously, and it is also possible to make the estimate with good computational efficiency.

Such an estimate may be generated in different ways. In particular, different formulas may be used for estimating different line properties.

Moreover, even though the estimate is made in dependence of Γ, it is not necessary to generate Γ explicitly. A related value may be generated and used, or the end result may be generated directly, without actually generating Γ as intermediate result.

For example, Γ, Γ² or one of them multiplied by a predetermined value, such as

$\frac{\Gamma^{2}}{\omega \cdot \hat{C}},$

or the real or imaginary part of any of those may be generated as an intermediate value.

According to one aspect of the invention, Γ may be generated by numerical solution of the equation Γ coth(Γ)=Z_(in)·jωĈ. According to another aspect, Γ² may be generated by solving a biquadratic equation that can be deduced by replacing Γ coth(Γ) in the expression Γ coth(Γ)=Z_(in)·jωĈ with a truncated Taylor series expansion of Γ coth(Γ) in the variable Γ. This is a computationally efficient way of determining Γ².

Other truncated series expansions or ratios of series expansions which correspond to Γ coth(Γ) may also be used.

As an option, accuracy of an estimate may be further increased by making several estimates of Z_(in) at different frequencies f, generating a set of intermediate values that correspond to the estimates of Z_(in), least squares fitting a curve to the set of intermediate values and estimating a property of the transmission line in dependence of the coefficients of the fitted curve.

An object of the invention is thus to provide a method for estimating at least one transmission property of a telecommunications transmission line.

An advantage of the method is that improved accuracy can be achieved. A further advantage is that computationally efficient implementations are possible.

Yet an advantage is that the estimates of Ĉ and Z_(in) may be made in many different ways, for example by measurement through a line interface board.

Another advantage is that the method may be used to improve on the method of WO2005SE01619.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a flowchart of a method according to the invention.

FIGS. 2-12 show flowcharts of examples of how the last step of the flowchart of FIG. 1 may be carried out.

FIG. 13 shows a diagram of an overview of the flowcharts of the FIGS. 2-7.

FIG. 14 shows a block diagram of an example setup for practicing the method of the invention.

FIG. 15 shows a diagram of frequency as a function of cable length for the condition that the absolute value of the propagation constant times line length equals π, for various ETSI cable types.

FIG. 16 shows a flowchart of a method for increasing the accuracy of estimates when several estimates of Z_(in) at different frequencies f have been made. This is an example of how the last step of the flowchart of FIG. 1 may be carried out.

DETAILED DESCRIPTION

With reference to FIG. 14, telecommunication transmission lines 60 typically terminate in one end at a customer premises site 61 and in the other end at a central office setting 130. Here, the lines are typically connected to line boards 100. The line boards 100 are typically controlled by a processor 110 of a telecommunication exchange. It is of great interest to be able to measure properties of the lines 60 from the central office setting 130.

A transmission line can be described by per-unit-length parameters: a series resistance R, a series inductance L, a shunt conductance G, and a shunt capacitance C (e.g in the denominations Ω/km, H/km, S/km and F/km respectively). These are called primary parameters of the cable type.

In practice, the shunt conductance of the transmission line can be neglected (G=0).

The complex-valued propagation constant γ and characteristic impedance Z₀ are defined in terms of the primary parameters: R, L, G and C, i.e. γ=√{square root over ((R+jωL)·(G+jωC))}{square root over ((R+jωL)·(G+jωC))} and Z₀=√{square root over ((R+jωL)/(G+jωC))}{square root over ((R+jωL)/(G+jωC))} respectively, where ω=2πf, f is the frequency and j is the imaginary unit. For an open-ended transmission line of length d, the complex input impedance can be expressed as

Z _(in) =Z ₀ coth(γd)

The values {circumflex over (R)}, {circumflex over (L)}, Ĝ and Ĉ, where {circumflex over (R)}=R·d, {circumflex over (L)}=L·d, Ĝ=G·d, Ĉ=C·d, are called line constants. Likewise, we introduce Γ=γ·d.

A method for estimating the line constants is described in patent application WO2005SE01619. There is a need to improve on the accuracy of the estimates generated with the methods described in that patent application to make them more useful, in particular the inductance estimate.

While the method of WO2005SE01619 uses series expansion in terms of line constants and angular frequency, the present invention uses in one aspect a series expansion in terms of Γ², and in another aspect a numerical solution for Γ. Contrary to WO2005SE01619 we can get one independent solution Γ²=Γ(ω_(k))² for each angular frequency ω_(k), instead of estimates of line constants for all frequencies. This reduces numerical problems related to that solution.

As a consequence, the present invention improves on the accuracy that is attained by the method of WO2005SE01619, but it can also be used together with other methods.

The accuracy improvement of the present invention is achieved as described below:

Neglecting shunt conductance G, (G=0), the propagation constant γ can be expressed in terms of the characteristic impedance Z₀ as follows

γd=Z ₀ ·jω·Ĉ

Since input impedance can be expressed as Z_(in)=Z₀coth(γd) one obtains

Z _(in) ·jω·Ĉ=Z ₀ ·jω·Ĉ·coth(Z ₀ ·jω·Ĉ)=Γ·coth(Γ)  (1)

where Γ=Z₀·jω·Ĉ.

If Ĉ, f, and Z_(in) at the frequency f are known, then Γ may be determined by solving eq. (1)

In general, a solution for Γ may be found by solving numerically equations of the type

${\coth (\Gamma)} = {{\frac{\Gamma_{in}}{\Gamma}\mspace{14mu} {or}\mspace{14mu} {\tanh (\Gamma)}} = \frac{\Gamma}{\Gamma_{in}}}$ where Γ_(in) = Z_(in) ⋅ jω ⋅ Ĉ

A numerical solution can be found using vastly differing techniques like e.g. the Newton-Raphson method or the van Wijngaarden-Dekker-Brent method. Those methods work best if the initial value is close enough to the true value of the solution. When solving for several different frequencies, a good guess can be using the solution from the previous frequency point as a start value for the next frequency point.

The following simple iterative methods may also be used for a numerical solution.

$\begin{matrix} {\Gamma_{i + 1} = {{arc}\mspace{11mu} {\tanh \left( \frac{\Gamma_{i}}{\Gamma_{in}} \right)}}} & (a) \\ {\Gamma_{i + 1} = {{arc}\mspace{11mu} {\coth \left( \frac{\Gamma_{in}}{\Gamma_{i}} \right)}}} & (b) \\ {{\Gamma_{i + 1} = {\frac{1}{2}{\log \left( \frac{\Gamma_{in} + \Gamma_{i}}{\Gamma_{in} - \Gamma_{i}} \right)}}}{where}{\Gamma_{in} = {Z_{in} \cdot {j\omega} \cdot \hat{C}}}} & (c) \end{matrix}$

The numerical solutions seem to work best for small values of Γ, such as is usually the case for frequencies in the POTS band. The simplest and most reliable numerical method appears to be the classic Newton-Raphson method.

However, we can create a polynomial to approximate eq. (1), which is entirely sufficient to solve the problem. Taylor series expansion of Γ coth(Γ) in the variable Γ, and limited to the first three terms gives the following polynomial equation for each frequency point.

$\begin{matrix} {{Z_{in} \cdot {j\omega} \cdot \hat{C}} = {1 + \frac{\Gamma^{2}}{3} - \frac{\Gamma^{4}}{45}}} & {{eq}.\mspace{14mu} (2)} \end{matrix}$

The series expansion is valid when the absolute value of line propagation constant times the line length is less than π, that is, |γd|<π or |Γ|<π. For a given frequency, this limits the length of the transmission line for which a property is to be estimated. Alternatively, for a given line length, the frequency is limited.

FIG. 15 shows a diagram of frequency as a function of cable length for the condition that the absolute value of the propagation constant times line length equals π, for various ETSI cable types.

The biquadratic equation according to eq. (2) is easy to solve and provides sufficient accuracy. However, it is also possible to use more terms of the series expansion, i.e. a higher order polynomial. General methods for finding roots (zero-crossings) of such polynomials are described in the technical literature.

Yet another equivalent series expansion may be derived from eq. (1) after expressing it as (Z_(in)·jω·Ĉ)·sinh(Γ)=Γ·cosh(Γ).

Replacing sinh(Γ) and cosh(Γ) with their series expansions and truncating after the three first terms, one may obtain after elementary manipulations

$\begin{matrix} {{{\frac{\Gamma^{4}}{120} \cdot \left( {5 - \Gamma_{in}} \right)} + {\frac{\Gamma^{2}}{6} \cdot \left( {3 - \Gamma_{in}} \right)} + \left( {1 - \Gamma_{in}} \right)} = 0} & {{eq}.\mspace{14mu} \left( {2a} \right)} \end{matrix}$

In general, coefficients of a biquadratic or higher order equation may also be obtained by other means, e.g. using elementary least squares techniques. Instead of using series expansion for cotangent hyperbolic, sinus hyperbolic and cosines hyperbolic, one can consider the function Γ·coth(Γ) as a polynomial ratio

$\begin{matrix} {{G(\Gamma)} = {{\Gamma \cdot {\coth (\Gamma)}} = \frac{b_{0} + {b_{1}\Gamma^{2}} + {b_{2}\Gamma^{4}} + \ldots}{1 + {a_{1}\Gamma^{2}} + {a_{2}\Gamma^{4}} + \ldots}}} & {{eq}.\mspace{14mu} \left( {2b} \right)} \end{matrix}$

The optimal values of coefficients b₀ b₁ b₂ . . . and a₁ a₂ . . . can be found by minimizing least square deviation from the function Γ·coth(Γ),

$\begin{matrix} {{V_{2}\left( {a,b} \right)} = {\sum\limits_{k = K_{1}}^{K_{2}}\; \left( {{\Gamma_{k} \cdot {\coth \left( \Gamma_{k} \right)}} - \frac{b_{0} + {b_{1}\Gamma_{k}^{2}} + {b_{2}\Gamma_{k}^{4}} + \ldots}{1 + {a_{1}\Gamma_{k}^{2}} + {a_{2}\Gamma_{k}^{4}} + \ldots}} \right)^{2}}} & {{eq}.\mspace{14mu} \left( {2c} \right)} \end{matrix}$

One possible simplification may be the assumption that a₁a₂ . . . are the coefficients of series expansion of the sinus hyperbolic; we can then find the coefficients b₀ b₁ b₂ . . . which minimize the least squares criterion

$\begin{matrix} {{V_{2}\left( {a,{b(a)}} \right)} = {\sum\limits_{k = K_{1}}^{K_{2}}\; \left( {{\Gamma_{k} \cdot {\coth \left( \Gamma_{k} \right)}} - \frac{b_{0} + {b_{1}\Gamma_{k}^{2}} + {b_{2}\Gamma_{k}^{4}} + \ldots}{1 + {a_{1}\Gamma_{k}^{2}} + {a_{2}\Gamma_{k}^{4}} + \ldots}} \right)^{2}}} & {{eq}.\mspace{14mu} \left( {2d} \right)} \end{matrix}$

Limiting our example to just a biquadratic equation, the problem is reduced to determining b₀ b₁ b₂ while

$a_{1} = {{\frac{1}{6}\mspace{14mu} {and}\mspace{14mu} a_{2}} = {\frac{1}{120}.}}$

A set of values of Γ └Γ_(K) ₁ . . . Γ_(K) ₂ ┘ covering an interval of interest is chosen, for example real values from 0 to 1 in steps of 0, 02.

For an improved match to Γ·coth(Γ), interval [K₁ . . . K₂] can be subdivided into several Γ intervals, where the estimation of b₀ b₁ b₂ will be performed. In practice, these intervals may correspond to values of Γ for e.g. low, medium and high frequencies.

Equations 2a-d and eq. 2 are some examples of how a biquadratic equation can be obtained.

A solution of eq. (2) will now be illustrated.

Eq. (2) has the form of a quadratic function of a square, as is shown by rewriting it

$\begin{matrix} {{{{- \frac{1}{45}}\left( \Gamma^{2} \right)^{2}} + {\frac{1}{3}\left( \Gamma^{2} \right)} + \left( {1 - {{Z_{m} \cdot j}\; {\omega \cdot \hat{C}}}} \right)} = 0} & {{eq}.\mspace{14mu} (3)} \end{matrix}$

Therefore the solutions of eq. (3) can be found using the classic quadratic formula:

$\begin{matrix} {\Gamma^{2} = \frac{{- \frac{1}{3}} \pm \sqrt{\left( \frac{1}{3} \right)^{2} - {4 \cdot \left( {- \frac{1}{45}} \right) \cdot \left( {1 - {Z_{m} \cdot {j\omega} \cdot \hat{C}}} \right)}}}{2 \cdot \left( {- \frac{1}{45}} \right)}} & {{eq}.\mspace{14mu} (4)} \end{matrix}$

This provides a solution with particularly good computational efficiency.

Observe also that

Γ²=(γd)²=−ω² ·{circumflex over (L)}·Ĉ+jω·{circumflex over (R)}·Ĉ  eq. (5)

This means that the accepted solution should have negative real part and positive imaginary part. The other solution is rejected.

From eq. (5), it can be deduced that an inductance estimate can be generated as

{circumflex over (L)}=−Re(Γ²)/ω² Ĉ  eq. (6)

And that a resistance estimate can be generated as

{circumflex over (R)}=Im(Γ²)/ωĈ  eq. (7)

In eq. (1), Γ was defined as Γ=Z₀·jω·Ĉ.

Hence, an estimate of the characteristic impedance may be generated as

$\begin{matrix} {Z_{0} = {\frac{\Gamma}{{j\omega} \cdot \hat{C}}.}} & {{eq}.\mspace{14mu} (8)} \end{matrix}$

Γ may be generated from Γ² by taking the square root of Γ² and selecting the root which has a positive real part and a positive imaginary part.

As shown above, then, when an estimate of Z_(in) at the frequency f, and an estimate of Ĉ have been made, an estimate of a transmission property of the line may be generated in dependence of Ĉ, f and Γ.

The methods disclosed may be combined in different ways to generate such an estimate.

For example, Γ may first be generated numerically. Γ is then squared to obtain Γ², and inductance and resistance estimates are generated by applying eqs. (6) and (7), (FIGS. 4 and 6) whereas an estimate of the characteristic impedance is made by applying eq. (8) to Γ (FIG. 2).

Alternatively, Γ² is generated by applying eq. (4). Eqs. (6) and (7) are then applied to generate inductance and resistance estimates (FIGS. 5 and 7), whereas for the estimate of the characteristic impedance, Γ is generated as the square root of Γ² having positive real and imaginary parts and eq. (8) is then applied (FIG. 3).

It is of course also possible to combine multiple equations into one when generating an estimate, as well as doing operations in many different, yet mathematically equivalent forms.

For example, one may note that the biquadratic equation (3) has a solution (eq. (5)) for which the real and imaginary parts have a common factor ω·Ĉ

Γ²=−ω² ·{circumflex over (L)}·Ĉ+jω·{circumflex over (R)}Ĉ=ω·Ĉ·(−ω·{circumflex over (L)}+j·{circumflex over (R)})

-   The common factor ω·Ĉ may be included earlier in calculations by     making a substitution

$\begin{matrix} {Y = \frac{\Gamma^{2}}{\omega \cdot \hat{C}}} & {{eq}.\mspace{14mu} (9)} \end{matrix}$

Equation (3) will be transformed into an ordinary quadratic equation

${Y^{2} - {\frac{15}{\omega \cdot \hat{C}} \cdot Y} + {45 \cdot \frac{\left( {{Z_{in} \cdot {j\omega} \cdot \hat{C}} - 1} \right)}{\left( {\omega \cdot \hat{C}} \right)^{2}}}} = 0$

If we denote by

$\begin{matrix} {p = {- \frac{15}{\omega \cdot \hat{C}}}} & {{eq}.\mspace{14mu} (9)} \end{matrix}$

and

$\begin{matrix} {q = {45 \cdot \frac{\left( {{Z_{in} \cdot {j\omega} \cdot \hat{C}} - 1} \right)}{\left( {\omega \cdot \hat{C}} \right)^{2}}}} & {{eq}.\mspace{14mu} (10)} \end{matrix}$

the solution is found by

$\begin{matrix} {Y_{1,2} = {{- \frac{p}{2}} \pm \sqrt{\left( \frac{p}{2} \right)^{2} - q}}} & {{eq}.\mspace{14mu} (11)} \end{matrix}$

Again we choose the solution with negative real part and positive imaginary part.

A resistance estimate is then generated as

{circumflex over (R)}=Im(Y)  eq. (12)

and an inductance estimate as

{circumflex over (L)}=−Re(Y)/ω  eq. (13)

An estimate of the characteristic impedance may be generated as

Z ₀ =√{square root over (−Y/ω·Ĉ)}  eq. (14)

The correct solution is the root having positive real part and negative imaginary part.

Alternatively, a resistance estimate may be generated directly as

$\begin{matrix} {\hat{R} = {{Im}\left( \sqrt{\left( \frac{p}{2} \right)^{2} - q} \right)}} & {{eq}.\mspace{14mu} (15)} \end{matrix}$

without generating Y. The root with a positive imaginary part is the correct solution, giving a positive value for {circumflex over (R)}.

An inductance estimate may alternatively be generated directly as

$\begin{matrix} {\hat{L} = {{{\pm {{Re}\left( \sqrt{\left( \frac{p}{2} \right)^{2} - q} \right)}}1\; \omega} + \frac{p}{2\omega}}} & {{eq}.\mspace{14mu} (16)} \end{matrix}$

without generating Y. The root which gives a positive value for {circumflex over (L)} is the correct solution.

By applying the rules of mathematics, numerous equivalent other formulas may be obtained which all embody the relationships disclosed above, explicitly or implicitly, and hence generate an estimate in dependence of Ĉ, f and Γ.

Reduction of the Influence of Noise, Etc.

In order to improve accuracy, estimates of Z_(in) at several different frequencies may be used. This is particularly convenient when Ĉ and Z_(in) have been produced by the methods of WO2005SE01619 since Z_(in) is then normally measured at a number of frequencies.

A simple approach is to make several estimates of a transmission line property based on values of Z_(in) at different frequencies f, and then average the estimates.

An approach which may provide greater accuracy is to generate a set of intermediate values that correspond to the estimates of Z_(in), and then to least squares fit a curve to the set of intermediate values and estimate a property of the transmission line in dependence of the coefficients of the fitted curve.

In particular, the real or imaginary part of Γ² may be the intermediate value. Several estimates of Z_(in) for different frequencies f are made, and for each such estimate, Im(Γ²(ω)) and Re(Γ²(ω)) are generated.

A straight line g(ω)=aω is least squares fitted to the value pairs of Im(Γ²(ω)), ω; and the line resistance estimate generated as {circumflex over (R)}=a/Ĉ.

A quadratic function h(ω)=bω² is least squares fitted to the value pairs of Re(Γ²(ω)), ω; and the line inductance estimate generated as {circumflex over (L)}=−b/Ĉ.

In general, after values of Γ² has been generated for a number of frequencies, this set of values may undergo noise reduction by replacing it with the function Γ²(ω)=(bω²+jaω). Estimates of transmission properties are then made using the values of this function in the same way as they would have been made from values of Γ² that had not undergone noise reduction.

In situations when additive frequency-dependent noise is present, which is normally the case, the weighted least squares method is preferably used to improve the accuracy of parameter estimation. Preferably weights are chosen that are inversely proportional to the noise variance at the particular frequency.

Modification of the Method of WO2005SE01619

If used with the method of WO2005SE01619, Z_(in) and Ĉ are preferably estimated as described in that patent application, while {circumflex over (R)}, {circumflex over (L)} and Z₀ may be estimated with the method of the present invention, making use of the estimate of Ĉ and the estimate of Z_(in) at a number of frequencies made with the method of WO2005SE01619. The calculations for the estimates of {circumflex over (R)}, {circumflex over (L)} and Z₀ are preferably made with the same equipment as is used for the estimates of Z_(in) and Ĉ according to WO2005SE01619.

DETAILED DESCRIPTIONS OF THE FIGURES

FIG. 1.

In the step 101, an estimate of the capacitance of the transmission line, Ĉ, is generated. The estimated capacitance is the total capacitance of the transmission line. The methods described in patent application WO2005SE01619 or other methods known in the art may be used to generate the estimate.

In the step 200, one or more estimates of the complex open-ended impedance, Z_(in), of the transmission line are generated. If more than one estimate are made, these may be made at different frequencies. For example, in WO2005SE01619 and the patent applications that it refers to, it is described a method to make an estimate of Z_(in) at a number of frequencies. Other methods known in the art may of course be used.

In the step 300, an estimate of a transmission property of the transmission line is generated in dependence of Ĉ, the frequency f and a value Γ, where Γ satisfies the relationship Γ coth(Γ)=Z_(in)·jωĈ, and ω=2πf.

The step 300 may be carried out in different ways. Examples of how the step 300 may be carried out are given in the FIGS. 2-12 and 16, and in the associated descriptions below.

FIG. 2. Numerical generation of Γ, estimate of Z₀.

In the step 301, Γ is generated by numerical solution of Γ coth(Γ)=Z_(in)·jωĈ.

In the step 305, an estimate of Z₀ is generated from Γ, ω and Ĉ by applying eq. (8).

FIG. 3. Calculation of Γ², estimate of Z₀.

In the step 302, Γ² is generated from Z_(in), Ĉ and ω by applying eq. (4).

In the step 304, Γ is generated by taking the square root of Γ². The correct root is the one with a positive real part and a positive imaginary part.

In the step 305, an estimate of Z₀ is generated from Γ, ω and Ĉ by applying eq. (8).

FIG. 4. Numerical generation of Γ, estimate of {circumflex over (R)}.

In the step 301, Γ is generated by numerical solution of Γ coth(Γ)=Z_(in)·jωĈ.

In the step 303, Γ² is generated by squaring Γ.

In the step 306, an estimate of {circumflex over (R)} is generated from ω, Ĉ and Γ² by applying eq. (7).

FIG. 5. Calculation of Γ², estimate of {circumflex over (R)}.

In the step 302, Γ² is generated from Z_(in), Ĉ and ω by applying eq. (4).

In the step 306, an estimate of {circumflex over (R)} is generated from ω, Ĉ and Γ² by applying eq. (7).

FIG. 6. Numerical generation of Γ, estimate of {circumflex over (L)}.

In the step 301, Γ is generated by numerical solution of Γ coth(Γ)=Z_(in)·jωĈ.

In the step 303, Γ² is generated by squaring Γ.

In the step 307, an estimate of {circumflex over (L)} is generated from ω, Ĉ and Γ² by applying eq. (6).

FIG. 7. Calculation of Γ², estimate of {circumflex over (L)}.

In the step 302, Γ² is generated from Z_(in), Ĉ and ω by applying eq. (4).

In the step 307, an estimate of {circumflex over (L)} is generated from ω, Ĉ and Γ² by applying eq. (6).

In the FIGS. 8-12, eq. (4) is applied implicitly, without generating Γ². The results are in dependence of the value Γ, even though Γ or Γ² are not explicitly generated.

FIG. 8. Calculation of Y, estimate of Z₀.

In the step 310, p is generated from ω and Ĉ by applying eq. (9)

In the step 311, q is generated from ω, Ĉ and Z_(in) by applying eq. (10)

In the step 320, Y is generated from p and q by applying eq. (11). As evident from eq. (9a), Y is equal to Γ² multiplied by a predetermined value.

In the step 330, an estimate of Z₀ is generated from Y, ω and Ĉ by applying eq. (14).

FIG. 9. Calculation of Y, estimate of {circumflex over (R)}.

In the step 310, p is generated from ω and Ĉ by applying eq. (9)

In the step 311, q is generated from ω, Ĉ and Z_(in) by applying eq. (10)

In the step 320, Y is generated from p and q by applying eq. (11). As evident from eq. (9a), Y is equal to Γ² multiplied by a predetermined value.

In the step 331, an estimate of {circumflex over (R)} is generated as the imaginary part of Y by applying eq. (12). As evident from eq. (12), the imaginary part of Γ² multiplied by a predetermined value is used for the estimate.

FIG. 10. Calculation of Y, estimate of {circumflex over (L)}.

In the step 310, p is generated from ω and Ĉ by applying eq. (9)

In the step 311, q is generated from ω, Ĉ and Z_(in) by applying eq. (10)

In the step 320, Y is generated from p and q by applying eq. (11). As evident from eq. (9a), Y is equal to Γ² multiplied by a predetermined value.

In the step 332, an estimate of {circumflex over (L)} is generated from ω and the real part of Y by applying eq. (13). Hence, the real part of Γ² multiplied by a predetermined value is used for the estimate.

FIG. 11. Direct estimate of {circumflex over (R)}.

In the step 310, p is generated from ω and Ĉ by applying eq. (9)

In the step 311, q is generated from ω, Ĉ and Z_(in) by applying eq. (10)

In the step 340, an estimate of {circumflex over (R)} is generated from p and q by applying eq. (15).

FIG. 12. Direct estimate of {circumflex over (L)}.

In the step 310, p is generated from ω and Ĉ by applying eq. (9)

In the step 311, q is generated from ω, Ĉ and Z_(in) by applying eq. (10)

In the step 341, an estimate of {circumflex over (L)} is generated from ω, p and q by applying eq. (16).

FIG. 13.

This figure is an overview of the FIGS. 2-7.

The figure shows e.g. that Γ is generated by numerical solution in step 301. In step 302 instead Γ² is generated by applying eq. (4). In step 303, Γ² is generated by squaring Γ, whereas in step 304, Γ is generated by taking the square root of Γ². The correct root is the one having positive real part and positive imaginary part. In step 305, Z₀ is generated from Γ, in step 306 {circumflex over (R)} is generated from Γ² and in step 307 {circumflex over (L)} is generated from Γ².

Hence, an estimate of Z₀ may be made by executing the sequence of steps 301+305 (corresponding to FIG. 2) or the sequence 302+304+305 (corresponding to FIG. 3). An estimate of {circumflex over (R)} may be made by executing the sequence of steps 301+303+306 (corresponding to FIG. 4) or the sequence 302+306 (corresponding to FIG. 5). An estimate of {circumflex over (L)} may be made by executing the sequence of steps 301+303+307 (corresponding to FIG. 6) or the sequence 302+307 (corresponding to FIG. 7)

As is clear from the figure, when more than one transmission property is to be estimated, all steps need not be repeated. For example, estimates of {circumflex over (R)}, {circumflex over (L)} and Z₀ may be generated by executing once the step 302 and storing the value of Γ², then executing the steps 304 and 305 for an estimate of Z₀, the step 306 for {circumflex over (R)} and the step 307 for {circumflex over (L)}, using the stored value of Γ². Alternatively, the step 301 may be executed once and the value for Γ stored, then the step 305 executed for an estimate of Z₀ and the step 303 to generate Γ², which is also stored. The step 306 is then executed for {circumflex over (R)} and step 307 for {circumflex over (L)}, using the stored value for Γ².

FIG. 14. Example of physical setting

As already mentioned, telecommunication transmission lines typically terminate in one end at a customer premises site 61 and in the other end at a central office setting 130. Here, the lines are typically connected to line boards 100. The line boards 100 are typically controlled by a processor 110 of a telecommunication exchange. It is of great interest to be able to measure properties of the lines 60 from the central office setting 130.

Hence, this is a setting in which the methods according to the present invention may advantageously be used. The generation of estimates of Ĉ and Z_(in) may, for example, be made by use of the methods disclosed in WO2005SE01619.

The calculations to generate estimates according to the present invention may be carried out in the processor 110, or in a separate computer such as a workstation 120. The results may for example be presented on screen, sent over a network, stored or printed.

FIG. 15.

This figure shows a diagram of frequency as a function of cable length for the condition that the absolute value of the propagation constant times line length equals π, for various ETSI cable types. The condition is valid below and to the left of each curve.

FIG. 16. Noise reduction of Γ² when several estimates of Z_(in) have been made, estimation of transmission property.

In the step 401, values of Γ² are generated for each of a number of estimates of Zin which were made at different frequencies f.

The generation of Γ² may for example be done for each frequency by numerical solution of Γ coth(Γ)=Z_(in)·jωĈ and then squaring Γ to obtain Γ², or by applying eq. (4).

In the step 402, a quadratic function h(ω)=bω² is least squares fitted to the value pairs of Re(Γ²(ω)), ω; and a straight line g(ω)=aω is least squares fitted to the value pairs of Im(Γ²(ω)), ω.

In the step 403, a transmission property is estimated, using for Γ² the function Γ²(ω)=(bω²+jaω).

If R or L is to be estimated, an arbitrary ω is chosen, Γ²(ω) is calculated and eq. (7) or eq. (6) is applied.

If Z₀ is to be estimated, the frequency f for which an estimate is desired is determined, Γ²(ω) is calculated for the corresponding angular frequency ω, Γ is generated as the square root of Γ²(ω)), and eq. (8) is then applied. The correct root is the one which generates Γ with positive real part and positive imaginary part.

Of course, even if several transmission properties are to be estimated, steps 401 and 402 need only be executed once.

Results

Early results indicate that the methods according to the present invention achieve estimates of {circumflex over (R)} and {circumflex over (L)} that may be 100 times more accurate than those of WO2005SE01619, provided that the effects of noise are suitably handled. 

1. A method for determining an estimate of at least one transmission property of a telecommunication transmission line, the method comprising the steps of: generating an estimate of a complex impedance Z_(in) of the transmission line at a frequency f; generating an estimate of a capacitance, Ĉ, of the transmission line: and generating an estimate of the at least one transmission property of the line in dependence of Ĉ, the frequency f and a value Γ, where Γ satisfies the relationship Γ coth (Γ)=Z_(in)·jωĈ, and ω=2πf; wherein a value X is generated from the relationship Γ coth (Γ)=Z_(in)·jωĈ, where X is Γ, Γ² Γ²/ωĈ or one of them multiplied by a predetermined value, or the real, or imaginary part of any one of those, and where X is used to make the estimate of the at least one transmission property; and wherein X is an independent solution for the angular frequency ω.
 2. (canceled)
 3. A method according to claim 1 wherein Γ² is generated as $\frac{{- \frac{1}{3}} \pm \sqrt{\left( \frac{1}{3} \right)^{2} - {4 \cdot \left( {- \frac{1}{45}} \right) \cdot \left( {1 - {Z_{in} \cdot {j\omega} \cdot \hat{C}}} \right)}}}{2 \cdot \left( {- \frac{1}{45}} \right)}.$
 4. (canceled)
 5. A method according to claim 1 wherein Γ is generated by numerical solution of Γ coth (Γ)=Z_(in)·jωĈ.
 6. (canceled)
 7. A method according to claim 5, wherein the numerical solution method is the van Wijngaarden-Dekker-Brent method.
 8. A method according to claim 5, wherein the solution for Γ is obtained by iterating according to any one of the functions Γ_(i+1)=arc tanh(Γ_(i)/Γ_(in))  (1) Γ_(i+1)=arc coth(Γ_(in)/Γ_(i))  (2) $\begin{matrix} {\Gamma_{i + 1} = {\frac{1}{2}\log \; \left( \frac{\Gamma_{in} + \Gamma_{i}}{\Gamma_{in} - \Gamma_{i}} \right)}} & (3) \end{matrix}$ where Γ_(in)=Z_(in)·jω·Ĉ.
 9. (canceled)
 10. A method according to claim 1 wherein the estimated property is line inductance.
 11. A method according to claim 10 wherein the inductance estimate is generated as —Re(Γ²)/ω²Ĉ.
 12. A method according to claim 1 wherein the estimated property is line resistance.
 13. A method according to claim 12 wherein the line resistance estimate is generated as Im(Γ²)/ω²Ĉ.
 14. A method according to claim 1 wherein the estimated property is line characteristic impedance.
 15. A method according to claim 14 wherein the estimate of the characteristic impedance is generated as $\frac{\Gamma}{{j\omega} \cdot {\hat{C}}^{*}}$
 16. A method according to claim 1 wherein several estimates of the same property are made using estimates of Z_(in) made at different frequencies.
 17. (canceled)
 18. A method according to claim 12, comprising the steps of generating further estimates of Z_(in) for different frequencies f; for each frequency f, generating Im(Γ²(ω)), ω; least squares fitting a straight line g (ω)=aω to the value pairs of Im(Γ²(ω)), ω; and generating the line resistance estimate as a/Ĉ.
 19. A method according to claim 10 comprising the steps of generating further estimates of Z_(in) for different frequencies f; for each frequency f, generating Re(Γ²(ω)); least squares fitting a quadratic function h(ω)=bω to the value pairs of Re(Γ²(ω)), ω; and generating the line inductance estimate as −b/Ĉ.
 20. A method according to claim 1, comprising the steps of generating further estimates of Z_(in) for different frequencies f; for each frequency f, generating Re(Γ²(ω)) and Im(Γ²(ω)); least squares fitting a straight line g(ω)=aω² to the value pairs of Im(Γ²(ω)), ω; least squares fitting a quadratic function h(ω)=bω² to the value pairs of Re(Γ²(ω)), ω; and making an estimate of a transmission property of the transmission line, using for Γ² the function Γ²(ω)=(bω²+jaω).
 21. (canceled) 